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In our statement of work for this grant, the first item is 
"support and coordinate with research underway at the 
Aerostructures Directorate . This interim report supports the 
design of the model rotor and the comparative study of coupled beam 
theory and the finite element analysis performed earlier at the 
Aerostructures Directorate by Robert Hodges and Mark Nixon. 

Attention is focused upon two matters — (1) an examination 
of the small discrepancies between twist angle predictions under 

t 

pure torque and radial loading and (2) an assessment of 
nonclassical effects in bending behavior. 

Our primary objective is understanding, particularly with 
regard to cause-effect relationships. Understanding, together with 
the simple, affordable nature of the coupled beam analysis, 
provides a sound basis for design. 

STATIC APPLIED LOADING CASES 

The three load cases considered by Hodges and Nixon have been 
considered here. The first case is bending due to lift and blade 
weight, the second is pure torque and the third is axial loading 
due to centrifugal force. 

There is some inconsistency in the equations for the applied 
loading. In the present work, the coordinate X is taken from the 
blade root, which is radial station 5.23. 

Bending Due to Lift and Blade Weight 

The distributed loading is 

q z = 0.02222X - 0.0123 (Ibs/in) (1) 
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The rotor model cross section appears in Figure 1. The 
coupled beam analysis of this loading case appears in Attachment 1. 

Beam deflection results appear in Figure 2. Bernoul i -Eul er , 
the classical engineering beam theory, results are denoted by "BE." 
This model is overly stiff. Also presented are three shear 
deformation models, SD1, SD2 and SD3, and the finite element 
results. 

The shear deformation model SI is an approximation obtained 
by setting the coupling stiffness and C 36 to zero. This is the 

classical shear deformation model in the spirit of Timoshenko. 
Clearly it is overly stiff also. This direct transverse shear 
effect is small for a beam of this slenderness. 

The complete theory, which includes all coupling effects, is 
denoted SD3. It provides good agreement with the finite element 
results. 

The approximation denoted SD2 is obtained by neglecting 
completely the classical shear deformation effect accounted for in 
SD1 in favor of the coupling mechanism associated with C 9 ^ and C 3 g. 
This model, therefore, includes only deformations due to the 
transverse shear-bending coupling and the usual bending 
contribution. The magnitude of this new, unexplored form of 
elastic coupling is seen to be enormous by comparing SD2 and BE 
results. This is a finding of major importance in understanding 
the behavior. 

The SD2 or SD3 models are required in this application in 
order to get sufficiently accurate predictions. This clearly 
excludes the earlier classical type theory of Mansfield and Sobey 
from practical use. 
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Pure Torque 


Although there was generally good agreement for the torsion 
case in the Hodges-Nixon comparative study, the effect of 
torsion-related warping was not included. This effect has been 
included in the analysis presented in Attachment 2. 

The classical St. Venant torsion theory result (without 
warping) is compared to the complete beam theory (CBT) and the 
finite element results in Figure 3. The CBT results, which differ 
from the classical (CL) only by the warping effect, are in 
excellent agreement with the finite element analysis. Restrained 
warping* creates a boundary layer zone near the blade root that acts 
to stiffen the blade and reduce the angle of twist. 

Axial Loading Due to Centrifugal Force 

This case is of the utmost importance because extension-twist 
coupling is to be used to control blade stall. In the Hodges-Nixon 
comparative study, the classical St. Venant theory was utilized for 
the coupled beam analysis. The discrepancy between analytical 
predictions and the finite element analysis was the greatest for 
this case. Classical theory was too soft and it overestimated the 
twist angle, a condition that is not conservative in view of the 
stated purpose of the model demonstration. 

As in the pure torsion case, the neglect of torsion-related 
warping is the reason for the discrepancy between coupled beam 
theory and the finite element analysis. A complete analysis of 
this loading case is given in Attachments 3 and 4. Attachment 3 
contains the overall response analysis. The axial force 
distribution is 

N = 913.83 - 7.875X - 0.75287X 2 (lb.) (2) 
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This expression differs from that quoted in the Hodges-Nixon work 
due to use of different coordinates. 

The twist angle distribution appears in Figure 4. The use of 
CBT brings the beam theory results in very good agreement with the 
finite element analysis. The rate of twist distribution is given 
in Figure 5. Again, the agreement is very good. 

Attachment 4 contains an analysis of the strain distributions 
for this loading case. The strain distributions are given in 
Figures 6 and 7. The results indicate that structural damage would 
be likely to occur at radial station 10 (X = 5) rather than at the 
root as predicted by classical theory. 

WARPING ANALYSIS 

A complete analysis of the effects of torsion-related warping 
appears in Attachment 5. Also included is a description of a 
simple warping model that is based upon a rectangular approximation 
for the cell. The equivalent rectangle is chosen to possess the 
same enclosed area. An assessment of this model suggests that it 
is adequate for the complete analysis. 

The main difficulty in accounting for warping is determination 
of the warping function and the stiffness Cjj. Both are 
accomplished readily with the approximate rectangular model. 


CONCLUSIONS 

In structures designed for extension-twist coupling, a high 
degree of bending-shear coupling is present which drastically 
causes the structure to be more flexible in bending. The impact of 
this effect on system performance must be assessed. 
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RADIAL STATION, IN. 
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Torsion-related warping is significant enough to warrant its 
inclusion in the beam analysis. A simple rectangular approximation 
may be used, which avoids the complexities associated with warping 
function and warping stiffness determination for sections similar 
to the D spar. With warping accounted for, the coupled beam theory 
is extremely accurate and easy to use. 
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Co. ordinates ...of the. . 0. ...spar .. under taken 

is 6h ocon. 


In. .. Fig. 1. in order to calculate the -carping 

function , 


one of the CO.ordina te , say ...JZ ., . JoaS to .be 

(rxjnd the 


junction of .... the. 

other one , g j la do so , 

the 

curve 


fitting method may be used . This annlgsd 

is. 

summarised 

Fa. Appendix I . 

--- - - - ■ ■ 




/In. . other . 

q|uick and UttCtfui approach. 

to 

find 

the 

cuarping ji mcHoa 

is to consider the D spar 

as 

tuuO 


region... .•, one. of . 

them b . rechznguiafr the 0 

tbcr 

on<z. 

IS 


a part of a circle. ... . This. approach, is shocun in. Fit)- 2 . 
























4 


The radius of fhe circhz. consider cd 


IS found ft? 


as 


R 


H 2 . L 2 

2H 


To 





[he fofal area . <4 1 + A 2 , is 

T + /Ao = R. arcsi>T. k 
1 3 2 R 

and is 

>4, s . k 

2 

•Tergfore Ihz area aJan+ed , /4, , is 

A. = jk ^rc'iiVik - (^-Hkk 

1 2 R 2 

The area of fHe reel-angular region , A 2 , Is 



L'.-H 


Ihus } Ihz lt>k<L enclosed ar<LQ of fhe D> spar is 


be 

C3) 
U€>2*d : 


C4) 

05) 

(lo) 

(7) 

(£>a) 


- 2 (_ 4 ^ + < 4 ^ ) 




b 



ORIGINAL FAGS & 

OF POOR QUALITY 


or. 

An. - 2 [ L / H. + — arcs ink. _ CR-H).k ] 

2 12 2 

(&b; 


.Tke circumference; of ... ike... circular, part region. , c 2 , 

is 


Cy - R . arc sin. k 

Z R 

(9) 


Tke circumference nf fh^„ rectangular region , £ r , is 



. C 1 , H + L / 

Oo) 

1-Ken 

the foha( circurrjferencc . of the D> spar, C, is 



c= .. .2 (c 1 + c 2 ) 

Clio) 

or 

c= 2 (_ L^+H 4 iR. oircsirv k ) 

(tlb) 


.Tkerefare kk- , . Lei sau 

C 



& 2 J*. 

02) 

So 

c 

L'H + arcsin.-^ - 

0 + M + R orQ sin k. 

' R - 

03; 


N 0 . 0 O , -fie ouarpin^ fun.ch’ari can be retonf-fzn 



cf-Cs) - B-S. _ 2ou($) 

- (W) 

c ni 

Since f ke main, purpose is fa find, .the Cuarpmg 
, cjii icla is calcu(af«:d by using a -Line, integral 

Cififfnesb .. 



K-h-M' 2 d * 


c 77 - $ 


(15) 


of different 

ma ter i oils - } 

H 


f C 77 ■ i 

CO 2 

foi 4d ds - 


and.noui, let consider +he lines shocon m Rg2. are made 
jt 

LdH C/2 


l'+H 
O^s i W 
U<: s $ L%M 


Sgm efr 


wKere 


^ = 65 - (S> 

4^ = &3 - 2caJ 2 (6) 

fo J = B 5 - 2 60 3 (s) 


l'+h $ s « | 


(17) 


(18) 


(19) 


so ( UJZ. reduced the prabl errs cohick ‘6 4a -f'n<d 4he 
sectorial areas swept out for eacU. region. . 

Fbr Line 1 


so 


Hron. 


r rt = l 


2^(5) . 1 L 'As 

O 

= L 5 


015^4 


(2t?) 
(21 ) 
(22) 


for Line 2 


fo = « 


so 


2ou 2 (s) = L'U + 


Hds 


(23) 

(24) 


4h 


en 


?w ? Cs)= l 4 + Hs )-| ^ s ^ L'+W 


(2b; 


7 


ORIGINAL PAGE 28 

OF POOR QUALITY 

For . Line. -3 , cue. nave to jioal. .. -the . shaded .. artza t A' , 

as juncfTan .. £?( 3 . The .’jarm . ..q^ A is (cund b d 
subs-bnacdng i he area bzlou) ih a shadizd. area ’f&.nv ike 
area of fhc circular ^>ort ... One., obtains the A 0i3 

A's .1 [ Rs _ C?.H).R Sin 1 j C26) 

2 R J 

tkzrC-(orz ._ ..._ ._. - 

\. ... 2cj 3 c$) 2L'U .* ?& _.XR-H),.Rsin_| l' + H £ (27) 

Using eguaHoni._C / '7 ) and £22) gireS 44 > 

. 4^(5) = (6- L')5 . tt (2&) 

£t.juaKonS C^&) ~and (25) . giv'CLS vf'j , . 

vk(s) = ( B-.HU. ~ L'H H 4t3.£L / +M . (29) 

And,. eguafions ..(19) ..and (?.?). gives . A 3 ., 

(&-R)i_.4(R.H)Eiki -.2L'H htW.2 S 2c ( 30 ) 

3 2 

Cgudkbn Ok) can be . cur il-fen 06 

/ on t*' \ . . 

c 77= 2( k M i 1 t ^ I 2+ r 3 J (31) 

c.v’huz.ie. . . .... 

n 

1 . a l[(B_L')5] 2 d'5 b2a) 


O 



8 


0F p ooR ou$!r£ 


-H 


/x* 


L = (b-OVda 


,.2 

K. LSlk2 S 3 

3 


11 


C 32 b) 


or 


-id 



L'+H 


I ? = 


[CB-H)5 _ l/H 1 2 c(S 


H 

L'+U 


I 2 . j[( v B-H)V 2 (b-H )i/H 5 + L'V]d 5 
H 


i'+« 


CHil s 3 _ S 2 + l'Vs 

3 


( 33 ) 


(34<0 


(340 


So 


<and 



c /2 


(B_R)s + (J2- H ) 12. sin — - 2L'W 


cti Gfo a ) 


L %4 


c /2 


•L 


(B_R)V * C?-H)Vsin 2 A + 4L'V + 2(0)^-H)Ps5,n 


t -44 


_ 4 O-<0Rl/W5in2 _ 4(B_R)L / H5 J d5 (3kdJ) 

? 


vnj (-4 
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ORIGINAL PAGE IS 
OF POOR QUALITY 




5 3 + ^_fl) 2 R 2 /l5 _ £ sin 25 \ 

3 ^ 2 4 R / 


4 R 

2 (&_R)(r.*OR (j2 2 5 in A . s.Rtosl) 

R R. 


- 4 CR-EORL'U (_£cot,A ) __ 2 


ZJz 




L+H 


1 3 > i 


A 3 _0_' + h.) 3 

8 


+ C*-«) * 2 )1(!-l'_h;_£ 


4 


C. ^;,,2(l 


sik; _ _ SiK 

EL 


t 44 h' 2 (5.l'_m) 

+ 2(V-K)(R-H)K 


R‘/s(kiA_ 


2R 


R 


A RcosA- _/L' + H ;Pco-. tA 

2 2R R 


f Kco^ Sl _ 

, PcoS L + W ) 

_ 2(V-K)l / U [ 

K 212- 

R ' 

L 


4 


SO c*J<2. UnouJ C<zcjUoh'or) {^^) } C ucaKon t35) j and 

1-5 (^i2Cjuothon ( 2 > 1 ) ) . 

■■If Ike ^ 5por 10 . mode. out of -ih& £>otma maler >d{, 

Oe>) 


k M . k m n ‘ ,5 K 

~ *^1-1 ■*■ "l-i ~ , 1'1 


Ik 


<Zr7 


C?? -:-2K n C l-| + l 2 + l 3 ) 


(39) 
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The Eq uivalent Ttecton g ular Area. A p proach 


Let consider a rectangular cross Section, which has 
-the same enclosed area as the D Spar . 

Linz 2. 

Lin<z 1 



It Ts presumed -that the rectangular has the same ajidth. } H 
aS the b Spar j, -there jorz 


<A) - 4bH 

t> span 


(.40 a ) 


or 


b = 


C 4 e ) 


£■ spar 


(40b) 


4ft 


So , the circumference: o£ -the reel angular Ci 0=0 section) 


10 


10) = 4 ( b + tt) 

rec 


C44 


-then the ratio , B ; 

r C 


B / = JML 


b + tt 


( 42 ) 


Tke jor« o.f the Suarpmg function j } 15 
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ORIGINAL PAGE IS 

0,7 POOR QUALITY 




Ers . _ 2cj ($) 


C43) 



We 

00 ill 

jbllocu 

. the 

Same procedure 0 * 

cehat 

cue have 

done 

jor 

b... 

spar 

case: 

. .1 .. 

- Let . soy Line 

1 

and 

Line 2 

•are. . 

made 

oui 

of. 

d fetar/t met lezrfa i s 


Uoi ng 

symmetry 

<3 he 5. 






.M _ 

fe 

+H 

r 




R 



C -0 2 

1 

r C2) 

2 ) 


— 


r 

u 77 ' 

- ^ • 

■ 4 

! K-i'i * ^ R1 ^ — *- 

■J 

1 

W do l (.44 ) 

R2 






c 

) 

H 



<4? here 


-- v f »■ 

1 

.3's • ._ 

. 0 66 £ H 

(45) 

. 'f 

1 R2 

Bj5 2oj 2 Co) 

k£ .s £ b + w 

(46) 

For Line . 1 

. ; 





r n = b 


( 47 a ) 

- - ■ 

2aJ^t*S)= j . baii 


b) 



0 

... . 2oj^(t>) = bs. . .. 

. ..0..&.S £ it 

(40) 


For . L \nz 2 


r K = H 


-50 


2oj^(^X » bH ^ J Hds 

o 

2 uj p (s) ; bH 4 Ha H ^ H + fc> 


( 49 a') 
b; 

C 60 ) 


Uoin0 - Qguabons C<A 5 ) an cl .(-4g>) gives 
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'Y = (&'- b )s 0^56 M 

1 El v 


Lc^u abono (4&) anol (do) g'iv<ze> 


^ = (B'_H)5 _ bH W45 6 b + H 


Equation (44-} can be r<zajni-Kzn a^> 


•R . C-n C2) . 

C77 - d (. **>i J-1 + ^2 ) 


coh&e. 


J, -- JUbE b)^] oi6 


0 

-w 


J 1= j U-b/5 2 cb = s~ 



and 


b+U 


2 


J, -- J l (p>Jw )5 _ bH J ds 


4 


b+H 


J 2 = j I (J 3 / - H )^s 2 — 2 (BifOb Hs + b ? B ? ]cb 

H , b+H 

(J 3 -H) %» ^BJd)bHs ? + b 2 M 2 .S 

3 


H 


bo 


J.. i®b!L>‘ 


( H+b) _ H' 


_ (b--H)bfi 


(b+bOiH 2 


or 


(51) 


( 52 ) 


(53 ) 

(34 ) 


J 55 ) 


( 5 b) 


) 4 b 3 ff 

J 
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AN APPLICATION! 


I yx our case- , Uac par ante \ar~ii c*re 

L". 0.5265 J 4=0,1469 i L= 0.0210 , R.1.3&G1 [In] 

Ikco IWc ' enclosed area, A e , and -line circumference , c, of” 
b spar are 

\4 e = a 2776- [in 2 ] i C= 2-6345 [in! 

\ t ■ 
f here -fore 

E>- 0.2107 ClMl 

Far [ +20 , "7 0 , -t 20, ~70 , '70, ■'2o 1 Graphite /Epo*y 

= 0.3139 £0 

and 

C 77 = 0.06&45 £5 [ lb- in 4 ] 

Usinc) the equivalent rectangular circa approacU 


if = 0-1469 ; b * 0.4724 


[tMl 


and 


4- 0-277(2 [iw 2 ] 


C= 2-4772 [IN] 


•So 


B / » 0-224-1 [ In] 


and 


C 77 - 0-05904 £5 [lb- in 4 ] 


r*. , 

"Tke error produced bz tuj<zer> C-j-j and C-j-j it> rJ 10.3/7 

but , -for /instance- tevist bp de pec-bod due -fa applied -fvrcjwL ca&z- 

44u error produced. [5 <v O.B % - 


